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A commentary on 

Force-field functor theory: classical force- 
fields which reproduce equilibrium quan- 
tum distributions 

by Babbush, R., Parkhill, J., and Aspuru- 
Guzik, A. (2013). Front. Chem. 1:26. doi: 
10.3389/fchem.2013.00026 

This comment regards a recent paper 
by Babbush, Parkhill, and Aspuru- 
Guzik (henceforth BPA) (Babbush et al, 
2013). The subject is the equilibrium 
thermodynamics of a system of many 
quantum particles with Hamiltonian H = 
Ef=iP-/(2m) + V(q), where q={q l ) 
and the commutator of the coordinate- 
and momentum operators is [q;, pj\ = 
ihhjj. In BPA it is correctly observed (see, 
e.g., Hillery et al, 1984; Cuccoli et al, 
1995) that the coordinate distribution 
function at temperature T = |} _ 1 , namely 



r\(q) = ^(q\e- m \q) 
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can be used to define the effective poten- 
tial W(q) (the kinetic part of the partition 
function Z is omitted in BPA) in terms 
of which the exact quantum equilibrium 
average of any operator 0(q) takes the 
classical form of a configuration integral, 

(O®) = J dq(q\0(q)e-^\q) 
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It is convincingly proven in BPA that 
the mapping r\(q) -o- W(q) (i.e., the 
exponential function) is bijective, as 
well as V(q) -o- r\(q), and then it fol- 
lows that V(q) -o- W(q) is one-to-one. 
Furthermore, it is correctly pointed out 
that the Giachetti-Tognetti-Feynman- 
Kleinert (Giachetti and Tognetti, 1985; 
Feynman and Kleinert, 1986) (GTFK) 
effective potential V e g(q) differs from 
W(q). Indeed, Veff(<?)> which accounts 
exactly for any quadratic potential, entails 
that for approximating (O(q)} one has 
to include a further Gaussian average 
accounting for purely-quantum fluctua- 
tions, as shown, e.g., in Vaia and Tognetti 
(1990); Cuccoli et al. (1995): there W(q) 
is also introduced and dubbed the local 
effective potential. However, in BPA it 
is not shown that Equation (29), the 
paper's main result obtained as the Jensen's 
approximation W(q) ~ W BPA (q) to the 
exact formula (26), can be calculated 
explicitly as 
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V[r(x)] 
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i.e., the convolution between the potential 
V(q) and a Gaussian with variance a 2 = 
fih 2 /(6m) proportional to the squared de- 
Broglie wavelength. This is in agreement 
with the Wigner series (Wigner, 1932) up 
to lowest order, but lacks the nonlinear 
contributions to it. How accurate is the 
approximation made in Equation (29) of 



BPA? One can estimate this by consider- 
ing a single (N = 1) quantum harmonic 
oscillator, V(q) = Kq 2 /2, whose frequency 
is ii) = y/vi/m. Expanding V (q + ^) in 
Equation (3) one finds 
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which does not improve upon the clas- 
sical result using Equation (2). From the 
known density for the quantum harmonic 
oscillator the exact functor for the class of 
harmonic potentials can be easily derived: 

moj 2 -j 
V(q) = —q 2 ^ 



1 , sinh fthtti 
W(q) = —In P 
2p pno) 



moo Bfcw 
+ — tanh— q 2 . (5) 



For a linear functor, this expression should 
be proportional to mm 2 : evidently this is 
true only in the classical limit, fihoo <C 1 or 
T 3> hw, where Equation (4) is recovered. 
However, the mapping V —> W can surely 
be locally linear, namely V + e&V W + 
e8W with 8W independent of the small 
parameter £. Hence, W SVA (q) is reliable 
only when the temperature overcomes 
the typical quantum energy scale ?ico; 
for instance, taking go 2 ~ V"(q m )/m [q m 
being the minimum of V(q)], a pair of 
hydrogen molecules has typically h(D ~ 
10 2 K (Vaia and Tognetti, 1990 and refer- 
ences therein) and using W BPA (q) would 
only be reliable at very high T» 10 2 K, 
i.e., just in the classical limit. Such an 
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approximation is indeed used in the high- 
T propagator of path-integral Monte Carlo 
algorithms in order to improve conver- 
gence in the Trotter number (Takahashi 
and Imada, 1984). Hence, the approxima- 
tion (29) of BPA can "reproduce quantum 
distributions" just when these are almost 
classical. 

On the other hand, the use of the exact 
effective pair-potential, rather than that 
obtained from Equation (29) of BPA, is 
a good starting point for treating a not 
too dense quantum fluid by means of a 
classical-like simulation, as shown in the 
last section of BPA and as noted by sev- 
eral authors (see, e.g., Thirumalai et al, 
1984 and many references cited in BPA). At 
variance with the procedure of BPA, based 
on the heavy calculation of a (locally) 
"linear functor" at fixed T, it would be 
more practical to directly obtain the exact 
pair-potential W(q) for the chosen V(q), 
a task that can easily be carried out at 
any T. 
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